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Abstract 

A general method for tight focusing of waves, based on compensation of the angular spectrum, is 
established. We apply the method to monochromatic, polychromatic and diffusive waves. Diffusive 
and monochromatic waves may form spatially localized waves in free space, whereas polychromatic 
waves form light sheets or non-decaying traveling evanescent modes confined to subwavelength 
regions in media where the frequency depends on the wave vector. We suggest an analogy between 
our compensation method and the transformation of frequencies between inertial, relativistic co- 
ordinate systems. 



I. INTRODUCTION 



Physicists and engineers often need to find methods for confining a field as much as 
possible with a minimum of knowledge of the system at hand. In the case of monochromatic 
waves, it has been found that phase and amplitude wave engineering is of crucial importance 
for subwavelength focusing of wavesfU E], |3]. In the case of pulsed waves, the situation is 
more complex, and optimal solutions have not been proposed to date, mainly due to the 
fact that the superposition of many spectral components in most cases tend to broaden 
the focal spot. Theoretically, one may hope to be able to focus strongly confined pulses as 
in Ref. m, but the aperture field and medium required remains elusive. However, there 
are a few promising pathways, e.g., using rainbow-colored apertures which combine to a 
single white-colored focus[5] or diffractive optical elements[6]. Well developed theories and 
numerical algorithms are available for studying monochromatic and polychromatic waves in 
focal regions in El 13 [9], but they do not provide direct insight into the reverse engineering 
problem and may not provide a fast solution or be optimal with the minimal amount of 
information given. Under such circumstances it is necessary to have straightforward methods 
which allows backtracking of the signal in a manner such that the initial field resulting in a 
strong focus is found. In this work we present a method to find the most strongly focused 
field. Here we extend the applicability of 'reverse' phase and amplitude engineering beyond 
previous approaches, and apply our method to monochromatic, polychromatic and diffusive 
waves. 

II. DIFFUSIVE WAVES 

We start by considering a related problem, namely the diffusion of a field in region of 
space containing no sources. The diffusion equation then describes the field according to 



where D is the diffusion coefficient. The field can now be written on the following general 
form 



where Ek{k) is the fourier amplitude. Usually, a solution of eq. [2] will represent a diffusion 
process where the field delocalizes with increasing time. However, we are here interested in 



dE{x,t) 



(1) 




(2) 



2 



finding a field E{x,0) with an optimally localized field E{x,tQ) {to > 0). To see how one 
may design such an initial field, let us assume that E{k) = EQexp[Dk'^to) is a constant for 
\k\ < ko and zero for \k\ > k^. Thus, the field is given by 

TP r^o 

E{x,t) = ^ e'''"-^'='(*-*o)dfc , (3) 
2vr J^ko 

The real value of eq. |4] at t = is given by 

TP rko 

E{x,0) = — cos{kx)e'''' '''dk , (4) 
27r J-ko 

whereas at t = to it is 

, ^ Eq [^^ sin{kox) , , 

E{x,to) = ^ / e'^^'dk = Eq , (5) 



2vr .Lko 



By selecting Dto such that Dk^to ^ 1, we see that the leading contribution to E{x,0) ~ 
cos{kQx)e^''o^° , which is an oscillatory function with period Aa; = 27r/fco- Thus, the field at 
t = is strongly delocalized, but with sharp local fluctuations, as seen in Fig. [l|a). On the 
other hand, E{x,to) is strongly localized with a central lobe 27r/fco, but is a factor e 



DkU 



smaller than the field at t = 0. A sharp localization is therefore obtained at the expense of 
an exponential reduction in field amplitude. 



III. MONOCHROMATIC WAVES 

The procedure introduced above can be straightforwardly applied to scalar monochro- 
matic waves satisfying the Helmholtz equation 

{V^ + k^) E{x,z) = , (6) 

where k = u/v is the wavenumber, u is the angular frequency and v is the velocity of light 
in the medium under consideration. Using the angular spectrum representation, the field 
can be written as 

E{x,z) = — / Ek{k,)e''^-+"'^'dk, , (7) 
where the angular spectrum is given by 

/oo 
0)e-^'^^^rfa;' . (8) 
-oo 
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where 



i^kl - F if P < kl 

Here k"^ > k'^ represent the homogenous plane waves, whereas k"^ < k^ correspond to inho- 
mogeneous or evanescent plane waves that propagate in the x direction and decay in the z 
- direction. 

The trick is now to express the angular spectrum as 

Ek{k,) = Eu{k^)e-'^^'^ , (9) 

where we must again require E{kx) to be bandlimited (i.e. nonzero only in a range ka < 
kx < kh) and converges sufficiently quickly such that the solution for the field does not 
possess any unphysical divergency problem. This was discussed in Ref. [11], and will not be 
repeated here. The field at the focal point {z = zq) is now given by 

1 

E{x,zo) = — / /fc(Me*'^"rffc, , (10) 

i.e. a Fourier transform of the aperture function /(fc^^) resulting in a finite resolution at the 
focal point. As an example, assume now that fk{kx) = 1 for \kx\ < ka {ka ^ k) and zero 



elsewhere. Then we have, as in Eq. 11 



E{x,Zoj = — / e dk = Eq ^ , (11) 
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In order to find the field E{x, 0) we must evaluate 
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E{x,0) = — I fk{kx)e"'^''-''''°dkx . (12) 

In Ref. [TT] we found that in the case of paraxial waves [k^ <^ k) the field E{x,0) is ap- 
proximately a quadratic phase function, which has been described in detail in standard 
textbooks on optics[l]. One can also design focusing of evanescent waves using this ap- 
proach. Let fk{kx) = Eq if ka < kx < kb and zero elsewhere. Here ka ^ k and Eq is a 
constant. The field can then be expressed as 

' ^iki,x-k,,(z-zo) _ ^ikaX-ka{z-zo)' 



E{x,z) ^ — 



(13) 



IX — {Z — Zq) 

The intensity at the focus is therefore just oc sinc^ ( ^''~^° x) whereas the aperture field 
oscillates rapidly within an envelope such that the intensity is oc l/(a;^ + z^). Thus, we have 



shown that the proposed method is in hne with the theoretical observations made in e.g. 
Refs. [31 [TT] for focusing of evanescent waves to subwavelength spots or hnes. Notice that the 
aperture field is much stronger than the field in the focal region, i.e. E{x, 0) / E{x, zq) ~ e*^*^". 
If we assume that ^ ka, a possible criterium for focusing could be obtained by requiring 
that the aperture intensity (at z = zq) envelope at x = zq, corresponding to the position 
with half the maximum intensity, should be wider than the focused intensity distribution 
of half- width ~ vr/A;;,. Thus, by requiring khZQ > tt, we ensure that the field is focused. At 
the same time, the field should not diverge at the aperture, and it is clear that kf,ZQ cannot 
be too large. Figure [l]b) shows the intensity / cx \E\'^ at z = Zq (solid line) and at z = 
(dashed line) when kb = I, ka = 0.1 and zo = 5. 



IV. POLYCHROMATIC WAVES 

An interesting question is whether the method above can be applied directly to polychro- 
matic waves (e.g. pulsed waves). To investigate this point, let us consider a time-dependent 
scalar field described by the wave equation as 

V^^(.,..O^i^!^ (14) 

In the angular spectrum representation, the field can be expressed as 

f POO POO ^ 

Fix, y,z,t) = Relj J E,it,, u)e'^^-+'^^-^-'^^-^^'^'dk^du | , (15) 
where k = u{kx)/c and 

^ ^ I v/^^ if k'>kl 

[ i^kl-k^ if P < kl 
In order to bring the polychromatic waves to focus we require that every spectral com- 
ponent is at focus a certain distance zi from the aperture at a given time ti. This can be 
done by compensating the time and space parts of the phase separately, or by requiring a 
cross-compensation. The first alternative requires that we set 

Ek{k,,u) = f{k,,uj)e-''^-^'^+'-'^ , (16) 

which results in 

E{x,y,z,t)= / /(A:^,tu)e*'=-^+*'=^(^-"^)"*"(*-*^)dMcu . (17) 



^-oo 



It is seen that the field in focus {zi,ti) is now just a sum over the frequencies of perfectly 
focused monochromatic wave components, which will be smeared out to a spot determined by 
the smallest wave vector of the spectrum (the largest wavelength in the case of homogeneous 
waves). In principle the procedure to find the aperture field is the same as that considered 
above for monochromatic waves. Unfortunately, it does not lend itself to immediate insight 
(except in the case of a rectangular spectrum for both kx and u, which results in similar 
fields as those given above), must in most cases be evaluated numerically, and will therefore 
not be pursued further here. 

The second alternative is to require cross-compensation of the phase, which amounts 
to setting kzZi — uj{kx)ti = 0, which basically states that the spatial part of the phase is 
compensated by the temporal part at {x,zi,ti). Such a requirement can only be fulfilled if 
u) depends on the spatial frequencies (i.e. direction of each plane wave), where the angular 
frequency is given by uj{k^) = i'y{vi)\kx\vi. Here 7(^1) = 1/ a/1 — (vi/cy and Vi = Zi/ti. 
The fact that frequency of the spatial wave vector depends on the direction was utilized in 
Ref. [3] to combine a rainbow spectrum to a single, focused spot. Conceptually the idea 
presented here is somewhat similar, although we use an entirely different approach to achieve 
the goal. Note that we must distinguish between the two cases vi> c and vi < c. 

Let us first consider the case vi > c. This can occur if the phase velocity of the wave 
exceeds the speed of light in vacuum. We then observe that 00 (k^) is real and given by 
^(kx) = akx {a = f 1/ a/ ("^i/c)^ — 1). The scalar field is then given by 



It is seen that the fields are time and space-shifted such that xj = {a/vi){z — Vit) gives 
transversal coordinate for the focal line. This suggests that the temporal field at {t = 
0,x,z = 0) can be imaged perfectly onto planes where z > 0. Our theory has implicitly 
assumed, as is the case for all real systems, a bandlimited angular spectrum, e.g. nonzero in 
ka ^ kx < kh. As an example, assume that Ek{kx-,uj) = Eq (i.e. constant) in ka < k^ < kb, 
such that the field is given by 



We notice that due to the peculiar dispersion relationship required by the phase matching, 
the waves propagate as a light sheet located at the plane z = {vi/a){x + at). At the 
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(18) 




(19) 



aperture, we require a pulse traveling in the negative x-direction. Oblique light sheets (in 
the case of electromagnetic waves) are then propagating outwards at an angle Vxja with 
the aperture. It should be noted that upon imposing the requirement ioik^^ = akj. without 
requiring additional boundary or media restrictions we obtain light sheets which somewhat 
unphysically extend to infinity (this problem could be removed by considering a more precise 
theory for causal and spatially dispersive media, but such an adventure is outside the scope 
of the current study). Only in the case Vi ^ c they are confined to the aperture plane at 
t = 0, but this requires that a small change in the spatial frequency generates a very large 
change in frequency, thus putting significant constraints on the medium to be use. 

A perhaps more interesting situation occurs when vi < c. Now a is purely imaginary, 
such that 

E{x,z,t)= [ [ Efc(A;^,Lj)e'''"^"'''"'^(*"^)dMw . (20) 

Jo J -oo 

As an example, we may approximate the angular spectrum as E{kx^uj) = Eq for ka < < 
kh 1 uJi < uj < UJ2 and zero elsewhere. Here we also assume that ka ^ fc, such that only 
evanescent waves are excited with kz ~ i\kx\ and oj{kx) ~ This evanescent-wave 



approximation therefore requires vi <^ c, and the field will be 

' ^iki,x—kt{z-vit) ^ikaX—ka{z—vit) 

E{x, z,t) oc 



ix — {z — vit) 



(21) 



Note that eq. 21 is identical to eq,13 if we set zq = vit, and is therefore just a translation 
of the evanescent wave along the optical axis with the intensity distribution similar to that 
of fig. [T] Since we require kb ^ k, which is fulfilled if we, as an example, set kf, = 10k = 
207r/A, the half-width of the intensity distribution at the traveling focal line is narrower 
than ~ 7c/kb = A/20. When t ^ l/k^v ~ ti our approximate theory for evanescent waves 
does not longer hold in the vicinity of z = (this follows from the considerations above; 
see also Refs. [IHl [TT] for a discussion about diverging solutions), where the intensity grows 
as exp{kbVit) / {x'^ + (fit)^) with time. The fact that the evanescent wave does not change 
at the focal line Zi = Viti is surprising, given the condition E{x,0,0) = E{x, Zi,ti) above. 
However, we also note that E{x,0,ti) ^ E{x, Zi,ti) such that the increasing energy of 
the field at the aperture is used to keep the field at the focal line zi = viti unchanged in 
magnitude as it propagates outwards. 

An interesting analogy occurs if one compares the problem of designing an aperture field 
that generate optimal focus with that of transformation of frequencies between inertial, 



relativistic coordinate systems[T2]. To see this, consider a plane wave of the form exp{kzZi — 
uj{kx)ti + (f){kx)). Any wavepacket is a weighted sum over such plane waves in spatial and 
frequency coordinates. In order to bring the wave packet to focus we require that every 
spectral component is at focus a certain distance zi from the aperture at a given time 
ti. This can be done by setting kzZi — uj{kx)ti + (pikx) = 0, such that the phase is exactly 
compensated at {zi, ti). It should be noted, as seen above, that such a requirement may give 
rise to strong localization at other positions as well. In any case, such a cross-compensation 
requirement can only be fulfilled if u depends on the spatial frequencies (i.e. direction of 
each plane wave), where the frequency is given by 



u{k:c) 



4>{k,)/t 



1 _^ 4>{kx)/ti 



1 + 



Vi\' 

c J 



- 1 



kxZi 
<P{kx) 



(22) 



Now consider the problem of transforming plane waves between inertial, relativistic coordi- 
nate systems. That is, consider an observer moving at a speed v relative to a fixed frame. 
According to Einstein's special theory of relativity, the moving observer will detect a fre- 
quency uj' given by uj' = 'y{v){uj — k^v), where 7(f) = 1/ \/l — (f /c)^P^. From the above 
it may be inferred that this expression corresponds to kzZ\ — uj{kx)ti + i4>{kx) = if we 
make the associations 4'{kx)/ti —* uj'/'~f{vi) and vi v, i.e. the phase factor (f){kx) must be 
associated with the frequency measured in the moving frame. In the special case considered 
in this study (pik^) = 0, which corresponds to u' = 0. For propagating waves we may set 
kz = {uj/cjcosO, where 6 is the real angle at which a specific plane wave makes with the 
direction of motion. Then we must have cos6 = c/v and therefore f > c, in agreement with 
the observations made above. 



V. CONCLUSION 



In conclusion, we have suggested a phase-compensation method for designing aperture 
fields giving rise to strongly confined waves. The idea is to first try to compensate the 
phase or amplitude such that only transverse spatial frequencies are left in the angular 
spectrum representation at a given region in space, thus allowing a strongly focused wave 
to form here. Next we calculate the aperture field required to obtain such a phase or 
amplitude compensation. The method has been employed to study diffusive, monochromatic 
and polychromatic waves, and shown to give new insight into the problems at hand. The 
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method here can probably also be applied to other wave systems where phase or amplitude 
compensation is beneficial. 
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FIG. 1: In (a) the diffusive waves (a) with ko = 1 (a.u.), D = 1 (a.u.) and to = 10 are displayed. 

In (b) the intensity of monochromatic, evanescent waves is displayed with ka = 0.1 (a.u.), kf, = 1 
(a.u.) and kbZo = 5 (a.u.). The solid lines show the field in focus, whereas the dashed lines show 
the fields at t = for diffusive waves and z = for monochromatic evanescent waves. 
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